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A finite-difference solution of the Schrodinger equation for negative donor centers D~ in two 

■ dimensions is presented. Our approach is of exact nature and allows us to resolve a discrepancy in 
the literature on the ground state of a negative donor. Detailed calculations of the energies for a 

, number of states show that for field strengths less than 7 = 0.1 17a. u. the donor possesses one bound 

■ state, for 0.117 < 7 < 1.68a.tt. there exist two bound states and for field strengths 7 > 1.68a.«. the 
Pj^ ' system possesses three bound states. Further relevant characteristics of negative donors in magnetic 

fields are provided. 
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The pro^^^paLjijlr centers in narrow quantum wells have become a subject of considerable interest during the 
past yearsaD'aDBD'ElBEI. In very narrow wells this system can be considered as a two-dimensional counterpartpQf 
the atomic H~ ion. A frequently used compound for current experimental investigations (see for example ref.c2l) 

S of such systems are layers of GaAs/AlGaAs. The high mobility i.e. small effective mass of the electrons and the 
comparatively large dielectricity constant of this semiconductor material allow us to study strong magnetic field 
effects in the laboratory. Both the two-dimensional character of the motion of the electrons and the external magnetic 
field makes the ground state of the system more tightly bound than that of the field-free threedimensional ion. 
Considering neutral donors D° the planar electron density of the 2d donor is enhanced significantly compared to that 
of the 3d donor. The form of the corresponding charge distribution makes the interaction of an additional distant 
charge with the neutral donor very different for the three-dimensional compared to the two-dimensional situation. In 
Q ■ three dimensions it is well-known that, for a fixed center (nucleus), the combination of a magnetic field witiuilie.|3d 
long-range attractive polarization forces lead to an infinite number of bound states for the negative donoiO-EJ'talii . 
In two dimensions the long-range interaction is of repulsive character. Indeed, it was shown by Larsen and McCanncI, 
I ' that an electron situated far from a 2d D° center experiences an overall repulsive potential which prevails both the 
^ . attractive Coulomb attraction due to the center and the attractive polarization forces. For distances of the order of 
7— I ' the extension of the neutral donor, exchange and correlation effects play a major role with respect to the binding 
ON , mechanism of the additional electron. 

Two-dimensional D~ centers in.qn agnjtic f j plds .. ar e described by a three-dimensional Schrodinger equation. The 
finite difference numerical methoc^OllZHlallj'E2lEllE3 allows us to solve this Schrodinger equation without any sim- 
plifications or approximations with respect to the geometry of the wave function, correlation effects etc. Due to 
this property of our computational approach it will be here possible to obtain a definite answer on the ground state 
energy and the number of bound states of the two-dimensional D~ center in magnetic fields. The critical values of 
the field strengths that are associated with the appearance of the bound states will be determined. One of the major 
motivations for the present work are the different theoretical values for the ground state energy ofp&-|two-dimensional 
1 ■ negative donor with and without magnetic field existing in the literature. Variational calculationsclu on the onehand 
and Monte Carlo simulations^] on the otherhand yield uncompatible predictions for the ground state energies. This 
puzzle will be resolved within the present work. 
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II. THE SCHRODINGER EQUATION AND THE METHOD OF SOLUTION 



The Hamiltonian of our two-dimensional system of two interacting electrons with an effective mass m and a singly 



charged positive ion placed into a magnetic field perpendicular to the plane can be written in Cartesian coordinates 
(x, y) as 



l,„o „ 7 s .if d d d d \ 7 2 

-j(Vf + V?) - .2 (-„_ +Xl —- y 2 — + x 2 _) 



±-L + ± (!) 



Here we use the units a e s = T^te. 2 m 1 = em* 1 • 5.29 x 10 9 cm for the length, E e s = e A h 2 me 2 — m*e 2 • 27.2eV 
for the energy and B e g = ce 3 m 2 h~ 3 e~ 2 = m* 2 e~ 2 -2.3505 x 10 5 T for the magnetic field strength. 7 — B/B c g, mo is the 
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mass of the free electron, m* = m/mo, e is the dielectricity constant of the semiconductor material, r\ — (x\ + y 2 ) 1 ^ 2 , 
r 2 = ( x 2 + an d r i2 = \ri — ^1- The spin terms are omitted. 

In contrast to variational calculations we are not biased by using a particular ansatz for the variational wave 
function but will perform a full grid solution to the Schrodinger equation that allows to control and estimate the 
remaining minor deviation from the exact eigenfunctions and eigenvalues. In the following we provide an outline 
of our approach which embodies various properties and transformations of the Hamiltonian and the corresponding 
Schrodinger equation. 

Transforming the electronic degrees of freedom to polar coordinates (fx, 4>\) and (r 2 , 4> 2 ) yields for the Hamiltonian 
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1 9 



1 9 

r 2 dr 2 



~5 'ax. 



r. 



1 9 2 



„-7 



' ( 2 



n r 2 



9 
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The eigenfunctions of this Hamiltonian are the eigenfunctions of the Hamiltonian (jl]) if we require them to be 
periodic with period 2ir with respect to the variables </>i and 4> 2 . The eigenfunctions of the Hamiltonian (0) should 
also be eigenfunctions of the z-projection of the total orbital angular momentum operator which is a conserved quantity 
ie. = Af* with 
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Introducing the coordinates ri,r 2 ,4>,& where 



and 
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we have 
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For these coordinates the Hamiltonian takes on the appearance 
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Eigenfunctions of the operator 1 2 read as follows 
We therefore have 
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This Hamiltonian depends on three degrees of freedom (n ,r%,4>) and the corresponding Schrodingpr.. equation can 
therefore be solved for the eigenfunctions and eigenvalues by applying our finite-difference methocLJ'E^lEJ. The coor- 
dinates r\ and r 2 are in the domain 

< ri < +00, < r 2 < +00 (11) 

with zero boundary conditions for infinity. The corresponding domain and boundary conditions for the angle <f> have 
to be analyzed in more detail. 

From equation (||) one can conclude that solutions of the Schrodinger equation belonging to the Hamiltonian jl(i| ) 
should be considered in the domain — 2ir < <fi < 2ir with cyclic boundary conditions. This leads to eigenfunctions 
that can be either symmetric or antisymmetric with respect to interparticle exchange. The first correspond to spin 
singlet states whereas the second yield spin triplet states. Beyond this we encounter additional eigenfunctions of 
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the operator (10) that are not eigenfunctions of our initial Hamiltonian (|l|) i.e. these eigenfunctions do not describe 
physical solutions. Indeed, exploiting certain symmetries of the wavefunction it is possible to confine the domain of 
the angle <fi to < <fi < 2ir for the Schrodinger equation belonging to the Hamiltonian ([l0|). As a consequence the 
non-physical solutions are excluded. For the states with M = solutions can be obtained in an even smaller domain 
< <f> < IT. 

The numerical finite-difference method employed here to solve the eigenvalue problem for t ha Ham iltonian (fLO) is a 
modification of the approach developed in previous works devoted to atoms in strong magnetidl3li3ll3 and electric iej 
(seeEJEil and references therein) and contains, in particular, technical aspects used in three-dimensional problemsEll 
Our computational procedure consists of the following main steps. The nodes have to be chosen in the spatial domain 
£1 : < n < +00, < r2 < +00, < <j> < 2ir (or < <j> < it). The values of the wavefunctions at the positions 
of the nodes are the numerical representation of the solutions of the initial differential equations. The nodes of the 
three-dimensional mesh in the space (n, r 2 , 4>) are placed at all the points with coordinates (ru, r 2 j, 4>k) where ru, r 2 j, 
and ipk are sets of the mesh node coordinates along the corresponding axes. For the coordinate 4> it is natural to use a 
uniform mesh with nodes at 4>k = 2ir(k — 1/2) /N for the region < <fi < 2ir (case M/0) and <fik = ^(k — 1/2) /N for 
the region < 4> < ir (case M = 0), where TV is the number of nodes in the direction 4>. For r\ and r 2 we have used 
non-uniform distributions of nodes similar to those described in ref.E3, which cover the infinite domains < 7*1 < +00, 
< r 2 < +00 as N — > 00. 

We employ the inverse iteration method to obtain the eigenfunctions and energy eigenvalues. This requires solving a 
system of linear equations with a matrix that is a finite-difference approximation to the Hamiltonian. The solution of 
these equations is particularly simple if the matrix has a block-tridiagonal structure. The latter can be achieved using 
the simplest three-point approximation for the derivatives for one of the coordinates. The choice of this coordinate is 
dictated by obtaining a convenient form of representation for the boundary conditions. For the two other coordinates 
we are free to apply more precise higher order approximations to their derivatives. The final values for the energy 
(and other observables) are provided by using the Richardson extrapolation technique for the corresponding results 
emerging from a series of geometrically similar meshes with different number of nodes. Using this approach we achieve 
a major increase of the numerical precision and, in particular, we obtain together with each numerical value a reliable 
estimate of its precisionEj'Lj. Typically, meshes used in the present calculations range from (the sparsest) mesh with 
20 3 nodes to (the thickest) one with 38 3 nodes, i.e. 38 nodes in each direction. 



III. RESULTS AND DISCUSSION 

Our results for the energies of the ground state of the two-dimensional D~ center are presented in table |. This 
table contains also the corresponding; energies for the ground state ofj-the D° center and a comparison with the best 
results obtained in the literatureaffil. RefBEl are variational and reffl represent Monte-Carlo simulations- One can 
see, that the variational results by Larsen and McCann coincide for weak and katermediate field strengths^ very well 
with our results. In contrast to this the results obtained by Louie and Tao Pangu differ significantly from these values. 
The absence of any approximations in our approach and the. possibility to reliably evaluate the convergence of our 
results allow us to conclude, that the results obtained in ref.Q overestimate the true values for the binding energy of 
the ground state of the two-dimensional negative donor considerably. (The reader should note that all digits of the 
values for our calculated energies given in table [j] are converged i.e. coincide with the exact results). This resolves the 
discrepancy on the ground state energy of the negative donor present in the literature as demonstrated by the results 
contained in table | particularly for 7 = but also for nonvanishing field strengths. 

For strong magnetic, fields 7 >> 1 table | provides also the energies obtained by Larsen and McCann in another 
variational calculationo. Our results both on the total and binding energies show that the ground state of the system 
is more strongly bound than predicted by ref.EI. The differences for the energies depend only weakly on the magnetic 
field strength and are approximately 1.08 a.u. for the total energy and approximately 0.36 a.u. for the binding energy! 
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TABLE I: Di 
[IS] and refsi 



nor total and binding energies in two dimensions for the M = singlet ground state. Results of our calculations 
[other]. The energies of the neutral donor are also given. 



7 


E D [IS] 


E D [other] 


E D [IS] 


E D [other] 


eW [is] 


eT\ 


3ther] 





-2.24027 


-2.239! 


-2.00000 


-2.0000 
-2.0000 


0.24027 


0.239E 
0.2555 






0.02 

0.05 

0.1 

0.2 

0.5 

1.0 

1.0 

2.0 

3.0 

4.0 

10. 

20. 

50. 

100. 

200. 

500. 

1000. 



-2.26014 
-2.28947 
-2.33712 
-2.42808 
-2.67354 
-3.02151 

-3.58733 
-4.05355 
-4.45883 
-6.27690 
-8.36994 
-12.5581 
-17.2937 
-23.9989 
-37.3107 
-52.316 



-3.0211 

-3.58ei 

-4.459;! 

-6.261? 

-8.3693 

-11.4583: 

-16.20452 

-22.916£^ 

-36.234; 

-51.243E 



-2.009967 
-2.024869 
-2.049518 
-2.098116 
-2.238416 
-2.455152 

-2.836203 

-3.165976 

-3.459582 

-4.815151 

-6.407114 

-9.62189 

-13.27207 

-18.449164 

-28.737120 

-40.33926 



-2.239 
-2.455 
-2.455 
-2.836 
-3.175 
-3.459 
-4.8133 
-6.405 
-8.8622 
-12.53314E 
-17.72454: 
-28.02495: 
-39.63327 



2tS 



0.25017 
0.26460 
0.28760 
0.32996 
0.43512 
0.56636 

0.75113 
0.88757 
0.99925 
1.46175 

I. 96283 
2.9362 
4.0216 
5.5497 
8.5736 

II. 977 



0.435 
0.566: 
0.585 
0.75. 
0.9 
l.OOOE 
1.463; 

I. 964E 
2.5961: 
3.6715: 
5.1923: 
8.2097: 

II. 61CE 



TABLE II: Total and binding energies for D triplet excited states with M < 0. 
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E M =-i 


E B 


- C/ M=-2 


E B 


E M =-3 


Eb 


0.05 


-2.02453 


-0.00034 


-2.02474 


-0.00013 


-2.02481 


-0.00006 


0.1 


-2.0492 


-0.0003 


-2.04916 


-0.00036 


-2.04933 


-0.00019 


0.2 


-2.100 


0.002 


-2.09715 


-0.00097 


-2.09759 


-0.000528 


0.5 


-2.253 


0.015 


-2.2360 


-0.0024 


-2.236554 


-0.001862 


1.0 


-2.50275 


0.04760 


-2.4523 


-0.0029 


-2.4511 


-0.0039 


2.0 


-2.9512 


0.1150 


-2.83775 


0.00155 


-2.8302 


-0.0060 


4.0 


-3.68474 


0.2249 


-3.472622 


0.01304 


-3.450805 


-0.008773 


10. 


-5.2634 


0.4482 


-4.85787 


0.04272 


-4.796 


-0.019 


20. 


-7.1028 


0.6957 


-6.4853 


0.0782 


-6.381 


-0.026 


50. 


-10.80161 


1.17972 


-9.7707 


0.1488 


-9.590 


-0.032 


100. 


-14.99309 


1.72102 


-13.501 


0.229 


-13.236 


-0.036 


200. 


-20.933 


2.484 


-18.794 


0.345 


-18.41 


-0.04 


500. 


-32.732 


3.995 


-29.25 


°' 51 J] 


-28.69 


-0.05 


500. 


-32.12890 


4.10090 


-28.62610 


0.6020 


-28.04640 


+0.0210 



The significant difference in deviations of our results from those obtained in ref.0 is due to the fact that the binding 
energies of the neutral donor in a strong magnetic field are underestimated in ref.H. As shown below this leads also 
to an overestimation of the binding properties of the excited states of the D~ center. 

The ground state of the D~ center is the only bound state for 7 = 0. We have carried out calculations for excited 
states of both symmetric and antisymmetric character of the spatial wavefunction with respect to the interchange of 
the coordinates of the electrons (i.e. spin singlet and triplet states). These calculations, were implemented both for 
M = and for M ^ 0. The latter states were investigated also by Larsen and McCannEl. Without loss of generality 
we confine ourselves to negative magnetic quantum numbers M < 0. Pairs of states being different only with respect 
to the sign of M possess equal binding energies within the corresponding Landau zones. However, opposite to states 
with negative M, states with M > do not belong to the lowest Landau zone. Therefore, it is reasonable to focus on 




values M < 0. 

Our calculations show that the ground state of the D~ center is the only spin-singlet state which is bound in the 
presence of magnetic fields. This conclusion coincides with results obtained in ref.B. On the otherhand, it follows from 
our calculations of the triplet states (see table ||) that there are two triplet states that become bound above some 
corresponding critical values for 7. These are the energetically lowest states for M = — 1 and M = —2, respectively. 
The M = -1 state becomes bound for 7 > 0.117 (the total energy of both D and D~ at 7 = 0.117 is -2.057852). 
The state with M = —2 becomes bound for 7 > 1.68 (the total energy at 7 = 1.68 is —2.720788). 

From table [n] it is evident that the lowest state with magnetic quantum number M — —3 is not bound for 
any magnetic field strength considered here and probably also for higher field strengths (the binding energy Eq = 
E — E D is negative and its absolute value increases with increasing field strength). This result differs from that 
obtained by Larsen and McCanro. They received a positive binding energies for this state for sufficiently strong 
magnetic fields. Their binding energy for 7 = 500 is presented in the last row of table, jnj. The main reason for this 
discrepancy is the underestimation of the binding energy of the neutral donor in ref.B as can be seen from table |. 
The data of table [0] together with our results for other excited states (they are unbound) allow us to conclude that 
the D~ considered possesses only two bound excited states - possessing magnetic quantum numbers M = — 1 and 
M = -2. 

Along with the total and binding energies of the two-dimensional negative donor D~ we have calculated some 
geometrical parameters of its wavefunction that provide additional information about the system. In figure [l] we 
present the expectation value R 2 — (r\ which characterize the extension of the spatial distribution of the 

electrons in the simplest and most straightforward way as a function of the field strength. First of all we observe 
that the electronic cloud is shrinking monotonically with increasing field strength for all bound states. This had to 
expected according to what we know about the behaviour of tightly bound states of few-electron system in strong 
magnetic fieldsEil. The difference of the behavior of R 2 for the ground (M = 0) and excited (M ^ 0) states for 
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relatively weak magnetic fields is obvious: The ground state is bound for 7 = 0, its wave function remains localized 
for all values of 7, and R 2 possesses a finite limit for 7 — ► 0. Furthermore it changes little for 7 < 0. la.it. For the 
excited states R 2 is not bounded for 7 — > 0. It changes rapidly with increasing field strength in particular in the 
weak magnetic field regime. Therefore R 2 possesses finite values for excited states only due to the presence of the 
magnetic field. On the otherhand, the dependencies of R 2 on the magnetic field strength in the high field regime are 
similar for all considered states since they are dominated by the diamagnetic term j(rf + r^)/^ of the Hamiltonian. 
The occurence of a small curvature in the dependence R 2 {^) for M = —1 between 7 = 0.1 and 7 = 10 is due to the 
influence of internal binding forces of the system, which are not negligibly small compared to the magnetic forces for 
this range of field strengths. 



IV. BRIEF SUMMARY 



We have shown that the two-dimensional negative donor D~ possesses three bound states ie. two bound excited 
states in the presence of a sufficiently strong magnetic field. The spin singlet ground state is bound for arbitrary 
field strengths. Our investigation of the ground state clarifies a discrepancy in.-lhe literature and confirms previous 
variational calculationsd'cl. As a result the energy eigenvalues obtained in ref.0 turn out to violate the variational 
principle and are too low i.e. the corresponding binding energies are too large. For weak and intermediate field 
strength we obtain good agreement of our total and binding energies compared to those of refJ3. In the high field 
regime, however, a significant lowering of the total energies and raising of the binding energies are obtained within 
the present investigation. Our particular computational approach allows for an estimate of the difference of our (up 
to several digits converged) results and the exact ones which consequently allows us to draw definite conclusions on 
the energies and properties of the donor. A series of calculations for excited states show that two other states become 
bound with increasing magnetic field strength. They are the lowest excited states with M = — 1 (7 > 0.117) and 
M = —2 (7 > 1.68). The extension of these states decreases monotonically with increasing field strength. 
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